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NON-EQUILIBRIUM CO CHEMISTRY IN THE SOLAR ATMOSPHERE 
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ABSTRACT 

Investigating the reliability of the assumption of instantaneous chemical equilibrium (ICE) for calcu- 
lating the CO number density in the solar atmosphere is of crucial importance for the resolution of the 
long-standing controversy over the existence of 'cool clouds' in the chromosphere, and for determining 
whether the cool gas owes its existence to CO radiative cooling or to a hydro dynamical process. Here we 
report the first results of such an investigation in which we have carried out time-dependent gas-phase 
chemistry calculations in radiation hydrodynamical simulations of solar chromosphcric dynamics. We 
show that while the ICE approximation turns out to be suitable for modeling the observed infrared 
CO lines at the solar disk center, it may substantially overestimate the 'heights of formation' of strong 
CO lines synthesized close to the edge of the solar disk, especially concerning vigorous dynamic cases 
resulting from relatively strong photospheric disturbances. This happens because during the cool phases 
of the hydrodynamical simulations the CO number density in the outer atmospheric regions is smaller 
than what is stipulated by the ICE approximation, resulting in decreased CO opacity in the solar chro- 
mosphere. As a result, the cool CO-bearing gas which produces the observed molecular lines must be 
located at atmospheric heights not greater than 700 km, approximately. We conclude that taking into 
account the non-equilibrium chemistry improves the agreement with the available on-disk and off-limb 
observations, but that the hydrodynamical simulation model has to be even cooler than anticipated by 
the ICE approximation, and this has to be the case at the 'new' (i.e. deeper) formation regions of the 
rovibrational CO lines. 

Subject headings: astrochemistry — molecular processes — radiative transfer — Sun: chromosphere 



1. INTRODUCTION 

Thirty years ago Noyes & Hall (1972) inferred very low 
brightness temperatures (T « 3700 K) from their discov- 
ery of strong rovibrational CO lines at 4.7 /im observed 
close to the edge of the solar disk. Several years later, 
after such surprising observational results had been con- 
firmed by Ayres & Testerman (1981), it was suggested that 
the low chromosphere might not be hot at all but could 
instead be permeated by CO-cooled 'clouds' at altitudes 
between 500 and 1000 kilometers above continuum optical 
depth unity (Ayres 1981). This led to controversy because 
other (UV and submillimeter) diagnostics had suggested 
the existence of a uniformly hot chromosphere with a min- 
imum temperature of about 4400 K near 500 km and a 
temperature rise above this temperature- minimum region. 
The controversy over the existence of cool gas in the low 
chromosphere continues today (see Kalkofen 2001, Ayres 
2002, Avrett 2003), after the publication of an abundance 
of literature on the subject with theoretical investigations 
concluding that the CO lines have LTE source functions 
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(Ayres & Wiedemann 1989; Uitenbroek 2000a) and with 
new observations showing off-limb CO emission protrud- 
ing hundreds of kilometers into the chromosphere (Solanki, 
Livingston & Ayres 1994; Clark et al. 1995), but also with 
the discovery that far-UV chromospheric lines observed on 
the solar disk always remain in emission at all positions 
and times (Carlsson, Judge & Wilhem 1997). 

Over the last few years, it has become increasingly ev- 
ident that the next crucial step towards a better under- 
standing of the enigmatic thermal structure of the so- 
lar chromosphere is to carefully investigate the reliabil- 
ity of the assumption of instantaneous chemical equilib- 
rium (ICE), which is currently used for calculating the 
molecular number densities in stellar atmospheres (see, in 
particular, Uitenbroek 2000a, b; see also Ayres & Rabin 
1996, Avrett et al. 1996, Asensio Ramos & Trujillo Bueno 
2003, Avrett 2003). Actually, both observations and sim- 
ulations indicate that the solar chromosphere is a highly 
inhomogeneous and dynamic region of low density plasma 
whose thermal, dynamic and magnetic properies we need 
to decipher for unlocking new discoveries. If the ICE ap- 
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proximation turns out to be adequate for modeling the 
strongest CO lines close to the edge of the solar disk, then 
the available CO observations would really be indicating 
the existence of cool gas in the solar chromosphere. Oth- 
erwise, a natural resolution of the current chromospheric 
temperature discrepancy could perhaps emerge if the CO 
concentration turns out to be sufficiently lower than would 
be expected on the basis of chemical equilibrium. This 
Letter addresses this challenging issue by showing the first 
results of a detailed investigation on the non-equilibrium 
CO chemistry in the solar atmosphere. 

2. FORMULATION OF THE PROBLEM 

Our strategy consists in performing chemical evolution 
calculations in the radiation hydrodynamical simulations 
of solar chromospheric dynamics described by Carlsson & 
Stein (1997, 2002), that do not include CO cooling in the 
energy equation. Therefore, at each time step of the hy- 
drodynamical simulation, we have fixed the ensuing ther- 
modynamic conditions and calculated the corresponding 
CO number density by following the chemical evolution 
starting from the molecular concentrations of the previous 
time step. 

The chemical evolution (CE) problem consists in calcu- 
lating the molecular number densities at each time step of 
the hydrodynamical simulation by taking into account all 
relevant chemical reactions with their ensuing rates. In or- 
der to solve this problem, one has to consider the following 
evolution equation for each chemical species T included in 
the model: 

it = Ea J2b Ec k ABcnAn B nc + Ea Eb k ABn A n B + 
+ J2a k A n A - J2a J2b kABinAnsni - 

-J2A k AinAn i ~k i n ij (1) 

where three-body reactions (first and fourth terms), two- 
body reactions (second and fifth terms) and one-body reac- 
tions (third and sixth terms) have been taken into account. 
When all the reaction rates for three-body (kABc), two- 
body {Uab) and one-body reactions (£u) are given, this set 
of equations for all the species included in the model form 
a set of first order ordinary differential equations which 
has to be solved numerically. Due to the stiffness of the 
system of equations (produced by the enormous differences 
in the abundances and their temporal variability), an im- 
plicit scheme has to be used. Two key ingredients have 
to be taken into account: the number of chemical species 
and the reaction rates for all possible reactions. After a 
careful investigation, we found that at least the follow- 
ing set of 13 species, which includes the most abundant 
diatomic molecules, is required to achieve a suitable de- 
scription of the relevant chemical processes: H, C, O, N, 
He, CH, CO, H 2 , OH, NH, N 2 , NO and CN. We have veri- 
fied that the inclusion of ionic species does not significantly 
affect the CO concentration for the typical atmospheric 
conditions encountered at heights /i<1000 km, although 
they should be ideally taken into account for very strong 
shocks capable of producing sizable changes in the degree 
of ionization. Concerning the reactions rates, we have used 
a reaction database created for the study of combustion 

2 This is similar to that used by Uitenbroek (2000a) and Ayrcs (2002) 
corresponding to the most dynamic cycle of the full simulation. 



mechanisms (Konnov 2000) which seems to be appropri- 
ate for the physical conditions in the solar atmospheric 
plasma. Other databases like UMIST (Bennet 1988) have 
been used for the study of the interstellar chemistry, but 
they suffer from the lack of certain important reactions 
which can efficiently take place in the Sun's atmosphere 
(e.g. the catalytic three-body formation of molecular hy- 
drogen: H+H+Hc^H 2 +He). We have also investigated 
the possible influence of CO photodissociations, which are 
one-body reactions. They are mainly produced by discrete 
photon absorptions at wavelengths between the Lyman 
cutoff (912 A) and the dissociation threshold (1120 A). 
By using the photodissociation rates of Mamon, Glassgold 
& Huggins (1988), we have verified that the contribution of 
photodissociation processes to the total CO concentration 
is negligible for the radiation field in the solar atmosphere. 

3. RESULTS AND DISCUSSION 

We used two time-series of snapshots from the above- 
mentioned radiation hydrodynamical simulations, each 
one lasting about 3600 seconds and showing the upward 
propagation of acoustic wave trains growing in amplitude 
with height until eventually producing shocks. The first 
one corresponds to a relatively strong photospheric distur- 
bance showing well-developed cool phases and pronounced 
hot zones at chromospheric heights (see Carlsson & Stein 
1997). The strongest dynamic cycle of this simulation pro- 
duces a peak-to-peak line-core brightness temperature of 
1000 K, concerning the strong 3-2 R14 CO-line synthe- 
sized at disk center with the ICE approximation. How- 
ever, the brightness temperature variations in most of the 
cycles of this simulation are 400 K approximately, which 
is similar to the observed values found by Uitenbroek, 
Rabin & Noyes (1994) under excellent seeing conditions, 
but larger than those inferred from the temperature his- 
tograms in Ayres & Rabin (1996). We shall refer to this 
as the strongly dynamic case. 2 The second simulation cor- 
responds to a much less intense photospheric disturbance 
(see Carlsson & Stein 2002), which for strong CO lines 
synthesized at disk center produces a peak-to-peak line- 
core brightness temperature fluctuation which is always 
smaller than 400 K. We shall refer to this as the weakly 
dynamic case. 

Starting from the molecular concentration given by the 
ICE approximation, we have followed the chemical evolu- 
tion in order to obtain the temporal variation of the CO 
number density (Nco) at each height in the simulated so- 
lar atmosphere. As expected, the CO concentration is an- 
ticorrelated with the local temperature variations, yielding 
relatively low Nqo values in the hot phases and relatively 
high Nco values in the cool phases. However, the am- 
plitude of the local Nco fluctuation is smaller than that 
given by the ICE approximation. In relatively low den- 
sity regions characteristic of the outer atmospheric layers 
(e.g. around h = 1000 km) the ICE approximation does 
a fairly good job during the hot phases, but it overesti- 
mates the CO number density during the cool phases. In 
contrast, in relatively high density regions characteristic 
of photospheric layers (e.g. around h — 400 km) the ICE 
approximation underestimates the Nqo values during the 

their ICE modeling, although they considered only a segment of 190 s 
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hot phases. Thus, the CO abundance does not react in- 
stantaneously to the changes in the temperature because 
of the finite reaction rates. If one carries out a linear anal- 
ysis by introducing a small temperature perturbation in a 
medium of given density and temperature and calculates 
the relaxation time needed to recover the original CO con- 
centration one finds that, for a given density, the relax- 
ation time is the larger the cooler the medium where the 
temperature perturbation is introduced. Similarly, for a 
given temperature of the unperturbed medium, the relax- 
ation time increases rapidly with decreasing density. Short 
relaxation times are typical of high-temperature and high- 
density media (e.g. t ro iax ~ 0.006 s for rin = 10 16 cm~ 3 
and T — 6000 K), while long relaxation times are charac- 
teristic of low-temperature and low-density situations (e.g. 
irda X «400 s for n H = 10 14 c m - 3 and T = 4000 K). One 
would find a broad range of relaxation times at a fixed 
height (e.g. h = 700 km) in the dynamic atmosphere that 
results from the above-mentioned hydrodynamical simula- 
tions, simply because the existing rarefactions, compres- 
sions and temperature fluctuations are continually chang- 
ing the atmospheric conditions. Obviously, the situation 
is highly non- linear and the relaxation time concept, al- 
though useful, loses its meaning in the real Sun. There- 
fore, any firm conclusion needs to be achieved via detailed 
numerical simulations. 

Fig. 1 refers to the strongly dynamic case. The left 
panel shows the temporal variation in the brightness tem- 
peratures of the line-core emergent intensities at \i = 1 
and fj, = 0.1 in the strong 3-2 R14 CO-line. The right 
panel gives the ensuing fluctuations of the atmospheric 
height where the line-core optical depth is unity, which 
we use as an indicator of the 'representative height' where 
the CO line-core radiation originates. As seen in the fig- 
ure, the ICE approximation does a sufficiently good job 
concerning the synthesis of the emergent CO spectrum at 
the solar disk center (fj, = 1), but it largely underesti- 
mates the line-core emergent intensities at fi = 0.1 dur- 
ing the cool phases of the hydrodynamical simulations, 
producing brightness temperatures that are typically 500 
K lower than those computed with the non-equilibrium 
CO concentration. Clearly, this is because during the cool 
phases the ICE approximation overestimates the 'heights 
of line formation' by about 300 km, concerning the synthe- 
sis of strong CO lines at /i = 0.1 in the strongly dynamic 
case. This happens because during such cool episodes 
the CO number density in the outer atmospheric regions 
is smaller than what is stipulated by the ICE approxi- 
mation, resulting in decreased CO opacity in the solar 
chromosphere. Interestingly, in the weakly dynamic case 
which has smaller kinetic temperature fluctuations (but 
still larger than the fluctuations of observed brightness 
temperatures!) the ICE approximation does a much better 
job even at /l* = 0.1. 

The reader may find surprising our conclusion that 
the ICE approximation is suitable for modeling the CO 
spectrum at the solar disk center, given that Uitenbroek 
(2000b) found that the spatially averaged line cores of 
weak CO lines synthesized in the three-dimensional (3D) 



hydrodynamical model of Stein & Nordlund (1989) are 
overly dark compared to the observed ATMOS spectrum 
described by Farmer & Norton (1989). As we shall show 
in a forthcoming publication (Asensio Ramos & Trujillo 
Bucno; in preparation), this is because such a 3D hydro- 
dynamical model of the solar photosphere is too cold in 
the CO line forming region. In fact, our ICE synthesis of 
the 7-6 R68 CO line in the improved 3D hydrodynamical 
model of Asplund et al. (2000) shows a notable agreement 
with the observed ATMOS spectrum, which constitutes 
an additional indication of the realism of the most recent 
3D hydrodynamical simulations of solar surface convection 
(see Shchukina & Trujillo Bucno 2001 concerning the iron 
spectrum in such a 3D hydrodynamical model). 

Fig. 2 contrasts the time-averaged CO concentration ob- 
tained assuming ICE at each time step of the strongly dy- 
namic simulation case with that resulting from the chem- 
ical evolution. Note that the ICE approximation leads to 
a significant overabundance of CO in the outer layers of 
the atmosphere (i.e. above 700 km). Thus, the CO num- 
ber density can be relatively low in such outer layers, in 
spite of the fact that the temporally averaged tempera- 
ture of the Carlsson & Stein (1997) simulations decreases 
outwards and does not show any chromospheric temper- 
ature rise. As expected from the previously reported re- 
sults, the ICE approximation does a sufficiently good job 
concerning the modeling of the temporally-averaged CO 
spectrum at the solar disk center (fj, = 1). In contrast, 
the emergent CO spectrum computed close to the edge of 
the solar disk (i.e. at fi = 0.1) shows significantly stronger 
CO-lines when the ICE approximation is used, especially 
concerning strong CO lines like the 3-2 R14 one, for which 
the line-core brightness temperature is about 100 K lower 
than that obtained using the non-equilibrium CO number 
densities. 

Finally, we turn our attention to the modeling issue of 
the off-limb CO emissions, which we have carried out by 
solving the radiative transfer equation in spherical coordi- 
nates at each time step of the hydrodynamical simulation 3 . 
The off-limb distances where (time-averaged) CO emission 
has been observed depend on the line (see Ayres 2002, for 
a summary of the available off- limb observational results): 
the off-limb emission extension of strong lines like the 3- 
2 R14 line lies between 0.55" and 0.7" above the 4.7 ^m 
continuum limb, while weaker lines like the 7-6 R68 line ex- 
tend to ^0.4". Fig. 3 refers to the strongly dynamic case. 
It shows that the atmospheric region where we can have 
a significant off-limb emission is extremely large when the 
ICE approximation is used, i.e. much larger than when the 
spectral synthesis is carried out using the non-equilibrium 
CO number densities. The dashed and solid lines show the 
corresponding height variation of the temporally-averaged 
off-limb emission at the core of the strong 3-2 R14 CO- 
line. They show that the non-equilibrium CO chemistry 
improves the agreement with the available off-limb obser- 
vations. First, as seen in Fig. 3, the ICE approxima- 
tion predicts that appreciable emission in the 3-2 R14 line 
should continue to relatively large off-limb distances, while 
partial eclipse measurements by Clark et al (1995) show a 



3 This is a suitable strategy for contrasting ICE and non-ICE results, as done also in Fig. 1. It is however clear that a truly realistic modeling 
of the available observations should take into account the geometrical foreshortening effect resulting from the fact that the solar chromosphere 
is a highly inhomogeneous medium in which the raypath goes through many different structures. 
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rapid disappearance of CO emission at heights above 700 
km, which is more in line with our chemical evolution cal- 
culations. Second, the representative off-limb emission ex- 
tension where the normalized intensity falls to 50% of the 
on-disk value is larger when the non-equilibrium CO con- 
centration is used, which helps to improve the agreement 
with the observations of translimb emission extensions, al- 
though they still seem to be below the most recently ob- 
served values (see Ayres 2002). 

4. CONCLUDING COMMENT 

Our results indicate that the CO line radiation observed 
close to the edge of the solar disk comes from atmospheric 
heights not greater than 700 km, approximately, and that 
the gas in these regions of the low chromosphere must be 
much cooler than indicated by the cool phases of the Carls- 



son & Stein (1997) hydrodynamical simulations. Lower 
temperatures will probably increase the relaxation times 
needed to reach the molecular equilibrium concentrations. 
In a forthcoming publication we will show what happens 
when the hydrodynamical simulations themselves arc car- 
ried out taking the CO cooling into account in a way con- 
sistent with the non-equilibrium evolution of the molecular 
number densities. 
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Fig. 1. — Temporal variation of the brightness temperature in the core of the strong 3-2 R14 line (left panel) and of the height of line-core 
optical depth unity (right panel) for the strongly dynamic case and for two observing angles: disk center (fi = 1, heavy lines) and close to the 
solar limb (fj, = 0.1, light lines), where fi = cosd (with 6 the angle between the solar radius vector and the line of sight). 
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Fig. 2. — Solid line: the height variation of the time-averaged CO concentration obtained from the chemical evolution calculation in the 
strongly dynamic simulation case. Dashed line: timc-avcraged CO concentration corresponding to the ICE approximation, but calculating 
the CO number densities of the atmospheric models associated to each time step by using the same chemical evolution code until reaching the 
ensuing equilibrium concentrations. Dotted line: time-averaged CO concentration corresponding to the ICE approximation, but calculating 
the CO concentrations directly from the Saha chemical equilibrium equations. A comparison of the dashed and dotted lines illustrates the 
reliability of the chosen database for the chemical evolution calculations. In any case, in order to be fully consistent with our comparisons, 
all ICE results in this paper refer to 'evolution until equilibrium' calculations. 
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Fig. 3. — The diagonal shading indicates the relatively large atmospheric region where the ICE approximation predicts off-limb emission 
in the core of the 3-2 R14 CO-line, while the horizontal shading shows that corresponding to the chemical evolution calculation. The figure 
also shows the calculated limb profiles for the continua at 4.7 (im and 5000 A . The zero level of the angular separation scale refers to the 4.7 
/im continuum edge calculated as the 50% point in the intensity drop, as normalized to the ensuing continuum intensity at disk center. At 
each angular distance, the vertical extensions of the shadings indicate the amplitude of the temporal variation of the line-core CO-emission. 
The dashed and solid lines give the corresponding temporally averaged off-limb emissions. The upper horizontal axis gives the height in the 
atmosphere, 2 = km being the location of continuum optical depth unity at 5000 A for a disk center observation (i.e., the usual depth 
scale in solar models). Note that at low altitude locations, where the IR continuum dashed-dotted curve lies above the (temporally-averaged) 
dashed and solid curves, the synthesized CO line turns into an absorption line for the ICE and non-equilibrium cases, respectively. 
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